Coherent Population Trapping with Controlled Interparticle Interactions 
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We investigate Coherent Population Trapping in a strongly interacting ultracold Rydberg gas. 
Despite the strong van der Waals interactions and interparticle correlations, we observe the persis- 
tence of a resonance with subnatural linewidth at the single-particle resonance frequency as we tune 
the interaction strength. This narrow resonance cannot be understood within a meanfield descrip- 
tion of the strong Rydberg-Rydberg interactions. Instead, a many-body density matrix approach, 
accounting for the dynamics of interparticle correlations, is shown to reproduce the observed spectral 
features. 
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Coherent population trapping (CPT), i.e. the pop- 
ulation of a quantum state decoupled from a resonant 
light field, serves as a paradigm for a quantum inter- 
ference effect First observed in 1976 @], CPT 
with its related phenomena electromagnetically induced 
transparency (EIT) 0, |1_ and stimulated Raman adia- 
batic passage (STIRAP) U has provided the basis for a 
large variety of effects and applications in many areas of 
physics, such as high-resolution spectroscopy, coherent 
control, metrology, quantum information and quantum 
gases. While CPT, EIT and STIRAP are generally de- 
scribed in terms of isolated single-atom interactions with 
coherent light fields, the situation becomes more involved 
when interactions between the particles need to be con- 
sidered. 

To gain initial insights into the effects of interactions 
on the quantum interference in CPT, consider two atoms 
with a three- level ladder structure with states |1), \2) and 
|3) as shown in Fig. Ufa). The atoms are exposed to two 
resonant coherent light fields and interact only if both of 
them are in the highly excited atomic state |3). In the 
case of non-interacting atoms the population accumulates 
in the two-body product state of the single-particle dark 
state \d) which is a coherent superposition of 1 1) and |3) 
[l[. This state is defined as the eigenstate of the to- 
tal Hamiltonian with vanishing coupling to the coherent 
light field. When turning on the interparticle interaction 
this state is no longer a dark state as it is no longer an 
eigenstate of the total Hamiltonian. As pointed out in 
[6j, the two interacting atoms, nevertheless, possess two 
dark states \d±). These states are dissipative due to the 
admixture of the intermediate, decaying state |2), but are 
significantly populated by optical pumping. While these 
states have dissipative character, they do not contain the 
state 1 33) and are, thus, immune to interactions. 

In a first approach to a many-particle system one could 
apply a meanfield model by replacing many-body opera- 
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FIG. 1: (a) Excitation scheme ( 87 Rb). fii and Q2 are the 
Rabi frequencies at 780 and 480 nm, respectively, 8 is the de- 
tuning of the upper transition; (b) calculated Rydberg state 
population, produced by the two-step sequence described in 
the text. The upper panel shows the result of a meanfield 
calculation, which predicts a strong shift and broadening of 
the resonance line. On the contrary, the exact result within a 
two-atom model (lower panel) yields an unshifted narrow in- 
teraction independent resonance and an additional resonance 
at S — V(R)/h with V(R) being the (repulsive) two-body 
interaction energy at an interparticle distance R. 



tors by products of their mean values thus neglecting in- 
terparticle correlations. Fig. [T^b) depicts the number of 
atoms in state |3) as a function of the upper laser detun- 
ing 5 and the pair distance, and compares the meanfield 
result (upper panel) to a solution of the fully correlated 
two-atom Bloch equations (lower panel). The meanfield 
model predicts a shift and significant broadening of the 
resonance line, due to the energy shift and decoherence 
of state 1 3) induced by the meanfield interaction with 
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the surrounding atoms. In contrast, the results of the 
two-atom model clearly indicate the population of the 
dark state as a non-shifted resonance without significant 
broadening. However, the states \d±) are non-separable 
and can therefore not be properly accounted for in the 
meanfield model. The blue-shifted resonance line in the 
two-atom model can be assigned to the weakly popu- 
lated eigenstate containing a contribution of 1 33) which 
is prone to interactions. 

In this Letter we address the question to which ex- 
tent a narrow unshifted CPT resonance persists in a 
multi-particle system with tunable interactions. More 
specifically, we investigate the three-level excitation of a 
gas of Rydberg atoms which are subjected to long-range 
van der Waals interactions in the blockade regime [7H9j. 
Rydberg atoms offer significant interparticle interactions 
over large distances which can be conveniently tuned by 
choosing the appropriate excited Rydberg state • The 
effect of the interactions can further be controlled by 
changing the density of ground-state atoms in the gas. 
Interaction-induced dephasing ll| and loss mechanisms 
[l2| on light propagation in a Rydberg gas under EIT 
conditions have recently been studied. The work pre- 
sented here provides a many-body approach which prop- 
erly accounts for interparticle correlations induced by the 
long-range Rydberg interactions. 

In our experimental scheme the states |1), |2) and |3) 
are represented by the atomic states 5Si/ 2 (i ? = 2, vnp — 
2), 5P 3 / 2 (F — 3,77j_f = 3) and a Rydberg state nS of 
87 Rb atoms which are coupled by resonant laser fields 
at 780 nm and 480 nm, respectively (see Fig.QJa)). The 
ground state atoms are trapped in a magneto-optical trap 
at a maximum density of po — 6.6 x 10 9 cm~ 3 . This 
maximum density of atoms in state |1) can be reduced in 
a controlled way by optically pumping atoms into another 
hyperfine state that does not couple to the excitation 
laser fields as described in Ref . [l3[ . The Rydberg atoms 
are detected by field ionisation. Further details on the 
experimental setup can be found in Ref. 

We apply a double pulse excitation scheme where the 
first pulse resonantly excites up to 20 % of the ground 
state atoms to the Rydberg state. The resulting well- 
defined mixture of atoms in the ground state and Ryd- 
berg state is probed by scanning the blue laser frequency 
of the second pulse. This initial partial excitation of the 
gas permits off-resonant excitation of strongly interact- 
ing atom pairs [lH during the second pulse, such that 
interaction effects are more pronounced. The first ex- 
citation pulse with a duration of 800 ns is realized by 
two circularly polarized counter-propagating laser beams, 
resonant with the respective transitions with peak Rabi 
frequencies of S7i = 7.6 MHz and O2 = 1.4 MHz, respec- 
tively. The lower (red) excitation beam has a large beam 
radius (wlmm), while the upper (blue) beam is focussed 
to a waist of «37 /im. The red Rabi frequency can thus 
be considered to be constant over the narrow cylindrical 
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ground state density p/Pbi 

FIG. 2: Excitation blockade of the Rydberg state 6fS (filled 
squares) in comparison to the excitation of the 30S state (filled 
circles) after the first (resonant) excitation pulse. While the 
30S state shows no saturation, the solid line in the 6f S data 
shows a heuristic saturation function tx 1/(1 + pbi/p) from 
which a blockade density pbi ~ 1-3 x 10 9 cm -3 is derived. 



excitation volume, while the blue Rabi frequency varies 
radially. After the first excitation pulse a second pair of 
laser pulses having the same beam geometry, but inde- 
pendently adjustable Rabi frequencies, probes the system 
with a pulse duration of 3 fis. While the lower laser tran- 
sition is still resonant (fii = 2.7 MHz), the upper probe 
transition is scanned over the atomic resonance with a 
Rabi frequency of VI2 = 1 .4 MHz. For the given parame- 
ters, on resonance the system is found close to the steady 
state after the probe pulse. A similar Rydberg excitation 
sequence has been employed in Ref. |15[ to probe energy 
shifts in a very dense sample with a detuned second exci- 
tation pulse. In contrast to our experiment, much shorter 
pulses were employed leading to a system far away from 
the steady-state. 

As a signature of interparticle interactions, the excita- 
tion blockade due to repulsive van der Waals interactions 
is presented in Fig. [5J There is no excitation blockade ob- 
served for the 30S state reflecting the n 11 dependence of 
the van der Waals interaction on the principal quantum 
number n [lfj |. 

The corresponding CPT spectrum for the 30S state 
is shown in Fig. [3ja). The observed CPT resonance 
can be well described in terms of single-particle Opti- 
cal Bloch Equations (OBEs) averaged over the Gaussian 
distribution of the Rabi frequency f^- The finite laser 
linewidth and redistribution of Rydberg states by black- 
body radiation have been included as additional decay 
processes. Consistency of measurements with the pre- 
diction of the OBEs was confirmed for various pulse se- 
quences and Rabi frequencies. Besides a scaling of the 
whole spectrum proportional to the density we do not 
find any density dependent features in the spectra. In 
the regime of strong Rydberg-Rydberg interactions (61S 
state with a van der Waals coefficient C§ = —1.18 x 10 21 
a.u. [n|), as shown in Fig. [2Jb) the CPT spectra still 
exhibit a pronounced peak with sub-natural linewidth at 
zero detuning at all densities. The width of the resonance 
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FIG. 3: (a) Probe scan of the 30S Rydberg state (upper graph) and simulation using one-atom Optical Bloch Equations (lower 
graph), averaged over the Gaussian distribution of 0,%. The 30S state is subject to stronger decay and redistribution than the 
61S state, which is why the dip in the signal around the central peak is less pronounced, (b) Upper graphs: Similar scans 
for the 61S Rydberg state at different densities. All densities are given in relation to the density pbi at which blockade effects 
become apparent. Middle graphs: Theoretical spectra, obtained from the density matrix expansion (red dots) and from the 
meanfield calculation (green dashed lines). Lower graphs: Theoretical spectra considering only interacting pairs. 



slightly increases with density as depicted by the squares 
in Fig. |U The CPT resonance width of » 3 MHz is well 
below the natural linewidth of the intermediate state \2) 
of 6.1MHz. As can be deduced from the comparison 
with the non-interacting 30S state (circles in Fig. HJ the 
resonance width is mainly determined by the finite laser 
linewidths. 

Our theoretical treatment of the excitation dynamics 
starts from the corresponding Heisenberg equations for 
the atomic population and transition operators a^l = 
\ a i)(fti\ ( a ,ft = 1)2,3), where i = 1,...,N labels the 
A^ atoms whose positions r 4 are randomly sampled from 
the underlying density distribution. Taking expectation 
values one obtains equations of motion for the reduced 
single-atom, two-atom, etc. density matrices, p^L — 



a/3/' Pa0,a'0' 



(frapfrrfpi) > respectively. The addi- 



tional interaction terms result in a hierarchy of coupled 
equations that ultimately requires knowledge of the N- 
atom density matrix for an exact solution and thus needs 
to be truncated in an appropriate way. The simplest 
possibility corresponds to the meanfield approximation, 

Pa0,a'0' - PafiPa'0' + 9 a 0,a'0' ~ Pa0Pa>0>> Le - Xne ne " 

gleet of direct two-particle correlations. Such meanfield 
treatments and slight modifications thereof 16 -IH have 
been successfully applied to model interaction-induced 
excitation suppression in cold Rydberg gases. While be- 
ing appealingly simple, the resulting nonlinear, single- 
atom equations, however, imply a unphysical distance- 
independent level shift for any finite Rydberg excitation, 
which has more dramatic consequences for the shape of 
the excitation spectrum. As depicted as dashed lines in 
Fig-EJb) the meanfield treatment, incorporating the den- 



sity profile as well as the spatial dependence of the Rabi 
frequencies of the experiment, fails once the excitation 
blockade sets in. 

To account for two-atom entanglement within a many- 
body description, the density matrix approach is ex- 
tended to second order, while approximately account- 
ing for three-atom correlations. This is achieved by 
truncating the hierarchy via the two-atom cluster ex- 

(Mj) (AO (W) i (*) ( fe j) 

pansion p\^p,^,p,, = P^Pa>p> , a » fi » + pJp>P a p,L"p" + 
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Pa"/3"Pa0,a>(S' *Pa0 Pa' 0' Pa" 0" ' 9 a0 ,a' 0' ,a" 0" 

tially neglecting direct three-atom correlations (i.e. 

9u0 l ,a'P' a"0")' one obtains a closed set of 37 dynamical 
equations for each pair of atoms. The dynamics of the 
two-atom density matrices involves three-atom interac- 
tion terms of the type 

E T/ (k,i,j) _ V"^ r T/ (i) (k,j) 

l jfe/'a l 8,a' i 3',a"/3" — [ V jkP a '0'Pa0,a"0" 



, (fc) , \r (j) (k,i) 

+ V ik Pa/3 9 a' 0', a" 0" + V jkP a " 0" 9 a 0,a' 0' 



(1) 



The second and third terms vanish rapidly if one of the 
three atoms is farther apart than the blockade radius 
from the remaining pair, but diverge if all three atoms 
are simultaneously very close to each other. This un- 
physical small-distance behavior results from neglecting 
direct three-atom correlations, which would cancel the 
respective terms. To account for this fact in a simple 
way we disregard the second and third interaction term 
in eq.fl}. This procedure, corresponding to the so-called 
ladder approximation in kinetic theory [19( , properly ac- 
counts for mutual interactions of close and distant parti- 
cles, and is approximate only if the mutual distances of 
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FIG. 4: The width of the narrow line in the measured 61S 
spectra (red squares) shows a weak density dependence and 
lies below the natural linewidth of the intermediate state 
at small densities. In comparison, the non-interacting 30S 
state does not show a significant change of the linewidth 
with density (blue circles). The many-body simulation (solid 
red curve) reproduces the slight increase of the linewidth for 
61S whereas the meanfield model (dashed green line) pre- 
dicts a much stronger density dependence. The black dash- 
dotted line shows the simulation results for vanishing laser 
linewidths. All values are FWHM determined from fitting 
the sum of two Gaussians of opposite sign. 



more than two atoms arc simultaneously on the order of 
the blockade radius. 

The solid dots in the middle graphs of Fig. [3Jb) show 
the results of this model qualitatively reproducing the 
experimental data. The calculation reproduces the weak 
density dependence of the resonance width (solid line 
Fig. [4j regardless of the strong excitation suppression ob- 
served for densities p > pu (see Fig. [5]). Again this stands 
in pronounced contrast to the strong quadratic broaden- 
ing obtained from the meanfield calculations. The sim- 
ulations also reveal that the measured resonance widths 
are largely limited by the spectral width of the excitation 
lasers as seen from comparison with the dash-dotted line 
in Fig. U giving the genuine CPT resonance width. 

The persistence of the narrow resonance may be qual- 
itatively understood from the two-atom dark state @ 
discussed above. In order to draw a connection to this 
simplified, yet more intuitive, picture and to elucidate 
the role of many-particle effects, we also consider the 
excitation dynamics for pairs of interacting atoms. As- 
suming that the interaction shift of each atom is solely 
determined by its nearest neighbor, one may calculate 
the pair-spectrum for variable distances, averaged over 
the nearest neighbor distribution of the random gas. As 
shown in Fig. |3] the pair model provides the basic mech- 
anism for the interaction-resistent narrow resonance. As 
may be anticipated from Fig. [1] the predicted linewidth 
becomes, however, density- independent for p > pbi, in 
contrast to the experimental observation. This is due to 
the oversimplified neglect of fluctuating, simultaneous in- 
teractions between several atoms, and more importantly 
interactions between blocked pairs, or larger numbers, of 



atoms. The improved agreement obtained using our re- 
duced density matrix approach, that accounts for these 
processes, thus, highlights the importance of many-body 
effects. 

In conclusion, we have presented an experimental 
scheme to investigate the role of interactions in a CPT 
scheme. As suggested by the occurrence of an entangled 
dark state in a two-body model, the interaction-induced 
correlations between the atoms cannot be accounted for 
by a meanfield model, and require a more sophisticated 
many-body theory. The theoretical framework developed 
in this work reproduces the observed density-dependent 
features of CPT resonances and complements current nu- 
merical methods that are limited to small samples and 
very small numbers of Rydberg excitations [17|, |2Cj, |21| . 
The concepts of CPT and EIT in systems exhibiting Ryd- 
berg interactions are currently attracting much interest 
in the context of quantum information, quantum simu- 
lation and sensitive probing of electric fields [f| 22 - 24 1 . 
Therefore the understanding of interaction effects is of 
crucial importance. The presented theory, supported by 
the reported experiments, may serve as a valuable basis 
for future studies of light propagation and EIT in inter- 
acting media, as well as addressing questions concerning 
optical nonlinearities and photon correlations induced by 
intcrparticle interactions. 

Note added- While this article was prepared for sub- 
mission we have become aware of complementary work 
[2o| , where the effect of interactions on EIT in a Rydberg 
gas is mapped onto the light field, leading to a coopera- 
tive optical non-linearity. 
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